Conditions for Adiabatic Spin Transport in Disordered Systems 
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We address the controversy concerning the necessary conditions for the observation of Berry phases 
in disordered mesoscopic conductors. For this purpose we calculate the spin-dependent conductance 
of disordered two-dimensional structures in the presence of inhomogeneous magnetic fields. Our 
numerical results show that for both, the overall conductance and quantum corrections, the relevant 
parameter defining adiabatic spin transport scales with the square root of the number of scattering 
events, in generalization of Stern's original proposal [1]. This could hinder a clear-cut experimental 
observation of Berry phase effects in diffusive metallic rings. 
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In contrast to phenomena related to Aharonov-Bohm 
(AB) phases [2] for charge carriers, the corresponding 
observation of Berry phases [3] due to the coupling of a 
spin to an orientationally nonuniform magnetic field B 
requires the limit of adiabatic spin evolution. In meso- 
scopic conductors, such a limit corresponds to the situa- 
tion where the carrier spin can follow the spatially vary- 
ing field during transport through the system. In terms of 
time scales, the adiabatic limit is reached when the Lar- 
mor frequency of spin precession, uj s = 2/iB/h, is large 
compared to the reciprocal of a characteristic time scale 
t c on which, from the point of view of the spin, the direc- 
tion of the field has changed significantly during motion. 
There is consensus that for ballistic (disorder-free) sys- 
tems with magnetic field configurations commonly the- 
oretically considered [1,4-6] and experimentally realized 
[7], it holds t c r~j L/vp, where v-p denotes the Fermi ve- 
locity of the carriers, and L is the characteristic length 
scale of the system over which the field changes. For 
one-dimensional (Id) ballistic systems the condition for 
adiabaticity, oj s ^> 2-K/t c , therefore reads [8] 
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where we introduced the adiabaticity parameter Qid- 

However, in the case of disordered systems there are 
two candidates for the characteristic time i c : (i) the mean 
elastic scattering time r and (ii) the Thoulcss time — 
(L/£) 2 t, with £ — vpT the elastic mean free path. The 
issue which of these two time scales is the relevant one 
has recently led to a controversial discussion [9,10]. 

In his proposal for Id diffusive rings Stern [1] perturba- 
tively calculated the lifetime of the adiabatic eigenstates 
and compared it to tiv He arrived at the condition 



Qid » L/e 



(2) 



for adiabatic spin transport (Eq. (7) in Ref. [1]). This 
corresponds to setting t c = r. Comparing Eqs. (1) and 
(2) one recognizes that in the diffusive regime, L ^> £, the 
adiabatic limit would require a magnetic field L/i times 



larger than in the ballistic case. This "pessimistic crite- 
rion", which later has also been advocated by van Lan- 
gen et al. [9] , would imply field strengths in the quantum 
Hall regime that let an experimental observation of Berry 
phases in diffusive metallic rings appear rather unlikely 
Alternatively, in analogy to the ballistic travelling time 
L/vp, it appears convincing to associate t c for diffusive 
systems with ixh, the time it takes the electron to dif- 
fuse through the structure. This argumentation has been 
put forward by Loss et al. [4,10,11]. By calculating the 
quantum corrections of the conductance in diffusive Id 
rings they predicted clear signatures of Berry's phase to 
be observable in a regime given by 



hd » e/L . 



(3) 



This condition for adiabaticity differs from criterion (2) 
by a factor (L/£) 2 and predicts for the observability of 
Berry phases a field strength above 20 mT [11], which is 
well in reach of modern experimental techniques [7] . 

In view of various recent experimental efforts to ob- 
serve Berry phases in the magneto conductance of meso- 
scopic rings [7,12-14], a clarification of this issue of the 
relevant time scale is desirable. The derivations of the 
conditions (2) and (3) were based on diagrammatic and 
semiclassical techniques. Here we choose a different ap- 
proach and study numerically the spin-dependent con- 
ductance of ballistic and disordered mesoscopic systems 
in the presence of a spatially varying magnetic field 
B(r) = V x A(r). The Hamiltonian for noninteracting 
electrons with effective mass m* and charge — e reads 
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The nontrivial coupling of the spin to the magnetic field 
enters via the Zeeman term fiB(r)-a, where a is the Pauli 
spin vector and fi = g*eh/ (Am^c) the magnetic moment 
with g* the gyromagnetic ratio. The electrostatic poten- 
tial V(f) includes the confinement and the potential of 
random impurities in the disordered case. At T = the 
spin-dependent conductance of a mesoscopic system with 
two attached leads is given by the Landauer formula [15] 
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Here i™ s m is the transmission amplitude from an incom- 
ing channel to with spin s to an outgoing channel (to', s'). 
We calculate <™ s m by projecting the corresponding Green 
function matrix onto the asymptotic spinors in the leads. 
We compute the Green function for the Hamiltonian (4) 
numerically, using a generalized version of the recursive 
Green function technique based on a tight-binding model 
[15] including spin [16]. We model the (non-magnetic) 
disorder potential leading to elastic scattering within an 
Anderson model by chosing random delta-like scatterers 
with amplitudes following a box distribution. The spin- 
dependent conductance is then obtained from ensemble 
averages over independent disorder configurations [17]. 

We now turn to the subject of interest and study how 
adiabaticity is approached in mesoscopic spin quantum 
transport. For this purpose we introduce a model system 
consisting of a 2d strip with a rotating in-planc magnetic 
field between two ballistic leads, see Fig. 1. This system 
can also be regarded as a model for transport through 
magnetic domain walls. We assume incoming electrons 
with spin-down polarization in the — y direction [18], in- 
jected from the left with Fermi wave number ftp = 2n/ Ap. 
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FIG. 1. 2d strip configuration used for the calculations of 
the spin-dependent conductance. The magnetic field B(x) 
performs a 180°-rotation within the plane of the strip. Spin 
states (Sin, Sout) are defined with respect to the y-axis. 

We first consider the overall conductance. In the bal- 
listic case, it is feasible to derive an analytical expression 
for the spin-resolved transmission of this system using a 
transfer matrix approach. The normalized transmission 
Tj^ for spin-down polarized incoming electrons to exit 
the system with spin-down polarization reads [19] 



1 



M 



M E 



\t m C-i\ 



M 

E 



m' ,m— 1 



M(i + Q 2 m ) 



(6) 



with the generalized adiabaticity parameter (Eq. (1)), 
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for the TO-th propagating mode in a 2d strip of length 
L and width W. Summing over all transverse modes 
in Eq. (6) we find that the overall dependence of the 



ballistic transmission (dashed lines in Fig. 2) is given by a 
Lorentzian ~ 1/(1 + Q 2 ) (dotted lines). This defines 
an 'effective' adiabaticity parameter Q ~ B for the 2d 
strip, with Qi < Q < Qm and Q <~ 1.4 Qid- This allows 
us to introduce a quantity that solely characterizes the 
adiabatic regime in the case of several open channels. 
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FIG. 2. Ensemble averaged normalized transmission (Tjj) 
for a disordered (solid lines) and ballistic (dashed) 2d strip 
(Fig. 1) as a function of the adiabaticity parameter Q ~ B. 
The panels correspond to different disorder strengths: (a) 
quasi-ballistic: L/t = 0.5 (L/W = 4.4), (b) moderate: 
L/l = 3 (L/W = 7.8), and (c) diffusive: L/l = 10 
(L/W = 13.9). The ballistic curves, Eq. (6), differ slightly 
from each other since they correspond to different Fermi wave 
vectors (ranging from krW/n = 7.7 to 11.6) but show the 
same overall Lorentzian decay with Q (dotted). 

For B — > the spin direction is preserved and Tj^ is 
maximal. In the limit of a strong £?-field the spin stays 
adiabatically aligned with the orientationally inhomoge- 
neous field during transport, minimizing the probability 
of leaving the conductor in Fig. 1 in a spin-down state. 
The Lorentzian dependence of on Q <~ B reflects 
this behavior and appears as the natural measure for the 
crossover from the non-adiabatic (Tj^ — > 1,Q <C 1) to 
the adiabatic (Tj^ -^0,Q>1) regime. 

To find a proper condition for adiabaticity in the disor- 
dered case we compute the ensemble averaged transmis- 
sion (Tn) in the presence of elastic scattering for Af <C £ 
as a function of Q and compare it to the ballistic re- 
sult (6). Our results for different ratios L/t are shown as 
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the solid lines in Fig. 2 exhibiting the following features: 
(i) The oscillations in the ballistic transmission are av- 
eraged out with increasing disorder, (ii) For (J » 1, 
the normalized (Tj^) is larger in the disordered than in 
the ballistic case. This means that in the presence of 
elastic scattering a stronger scaled field Q is required 
for acceding to the adiabatic regime of (Tj^) w 0. (iii) 
For Q < 1 we observe the opposite behavior: The non- 
adiabatic limit of almost maximum transmission (T^) 
is restricted to lower magnetic fields compared to the 
ballistic case, (iv) The crossover region Q ~ 1 is char- 
acterized by a transmission plateau, which approaches 
(I'll) ~ (^Ti) ~ 0.5 with increasing diffusiveness. Here, 
the non-magnetic disorder acts as a spin randomizer of 
the originally spin-polarized current. The features (i)- 
(iv) already begin to appear in the quasi-ballistic regime 
(Fig. 2(a)) and become more pronounced with increasing 
degree of diffusiveness given by the ratio L/£ (Fig. 2(c)). 

After this qualitative discussion we now derive a quan- 
titative condition for adiabaticity in the disordered strip. 
In analogy to the ballistic case we expect the disorder av- 
eraged transmission (Tn) to exhibit a scaled Lorentzian 
dependence in the limits of small and large Q. Indeed, 
in the adiabatic limit the Lorentz function is an excellent 
fit to the diffusive curve (T±±) in Fig. 2(c), if the ballis- 
tic parameter Q is replaced by Q/\fN a ^ with iV ac i being 
fitted (left inset Fig. 3). Corresponding results hold for 
the nonadiabatic limit where we use the scaling QVMiad 
(right inset in Fig. 3). We further determined iV a d and 
A^nad for various ratios L/£ and obtain power-law depen- 
dences illustrated in Fig. 3. We hence can formulate as a 
necessary condition for adiabatic spin transport through 
the disordered 2d strip: Q 3> (L/£) - 95 . Comparing this 
with Eq. (2) we obtain a little smaller exponent. To ex- 
plain this deviation we note that Eq. (2) can be written 
in the more general form Q » y/ (N), with the average 
number of scattering events (N) = (tTh/fr = (L/£) 2 . 
This suggests to associate _/V a d and iV na <j with the num- 
ber of scattering events the electron has to undergo upon 
traversing the microstructure. Due to the strong cou- 
pling of the finite-size 2d strip to the ballistic leads we 
expect the diffusion time to be smaller than the Thouless 
time frrh, thus reducing the number of scattering events. 

To confirm the above arguments we independently 
checked numerically the dependence of (N) on the scaled 
length L/£ of a finite disordered conducting strip. To 
this end we used a Id [20] random walk model taking 
explicitcly into account the interface between disordered 
and ballistic regions. We find that (N) obeys a power 
law in L/£ with exponent 1.91 which, as expected, is 
lower than two (dotted line in Fig. 3). In Fig. 3, besides 
small deviations for small Lji in the non-diffusive limit, 
there is good agreement with the fitted straight lines for 
Lji > 5, indicating diffusive behavior. Within the given 
error tolerance all three curves in Fig. 3 exhibit identical 
exponents and deviate only in the prefactor of order one. 



We conclude from our numerical, quantum mechanical 
results, together with the expression for (N) from the in- 
dependent random- walk model, that the adiabaticity pa- 
rameter scales with (N). This enables us to formulate a 
general (system-independent) adiabaticity condition for 
diffusive systems, that only depends on the correspond- 
ing adiabaticity parameter Q of the ballistic system and 
the mean number of scattering events (N) [21]: 

Q » ■ (8) 

For diffusive Id rings this criterion is in perfect agreement 
with Stern's original condition (2). 
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FIG. 3. Functional dependence of fit parameters iV a d 
(solid line) and iV na d (dashed) on the scaled length 
Lji of the disordered strip. Linear regression yields 
iV ad = 0.31 [Ljt) 1 - 9 (solid line) and 7V nad = 0.18 {L/t) l ss 
(dashed). For comparison the mean number of scattering 
events (N) = 0.48 (Lji) 1 ' 91 is also shown (dotted) obtained 
from an independent Id random walk model. Insets: Trans- 
mission {Tii) an d fitted Lorentzians in the nonadiabatic (left) 
and adiabatic (right) limit for strip with Lji = 10. (Error 
bars include uncertainty from fitting procedure.) 

So far we considered the total conductance which is 
dominated by the Boltzmann contribution. However, sig- 
natures of Berry phases in diffusive conductors appear 
only in the phase coherent part of the conductance, i.e. 
quantum corrections such as Aharonov-Bohm (AB) oscil- 
lations and universal conductance fluctuations (UCFs). 
To decide whether distinct Berry phase effects, e.g. in 
diffusive rings, can be observed at realistic magnetic field 
strengths, one has to check if an adiabaticity condition 
different from (8) holds for the quantum corrections. For 
this purpose, spin-resolved UCFs represent a suitable 
quantity, defined as 5g s / s = \J (T 2 , s ) — (T s / S ) 2 in units 
of e 2 /h. We calculated 5g s > s numerically as a function of 
Q for a diffusive 2d strip with L/£= 15. The results are 
depicted in Fig. 4 in terms of the normalized difference 
{^9ll~^9ll)/(^9ll+^9ll) which can be regarded as a po- 
larization. We note that the UCFs exhibit precisely the 
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same scaling behavior as the corresponding quantity for 
the total conductance, (Tn — 7ri)/(^n + ^Tl)> an( ^ con " 
sequently obey condition (8). Fig. 4 also illustrates the 
important fact that in a wide region £/L < Q < Ljl the 
respective adiabatic (T^, 5g^[) and nonadiabatic (Tj^, 
^5u) components are comparable in magnitude. 
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FIG. 4. Quantum fluctuations 5gn/(5gn + Sgn) (dashed) 
and Sg^ i/(Sgn + 5gn) (dotted) as a function of Q for a dif- 
fusive strip with L/l — 15. For comparison we also show 
{$9ll ~^9n)/(^9ll +$9u) (dashed-dotted) and the polariza- 
tion «T U > - <T n »/«T u ) + (T Ti )) (solid). 

We further note that numerical quantum calculations 
of the spin-dependent magneto-conductance in disor- 
dered rings subject to a circular inhomogeneous i?-field 
show that signatures of Berry phases appear only in the 
AB-oscillations of the adiabatic components ((Tri) an d 
Sgti), which are dominated by electrons with spin always 
aligned with the local field [22]. In view of Fig. 4 we 
hence can conclude for diffusive rings that in the exper- 
imentally relevant plateau region £/L < Q < 1 the adi- 
abatic components, which show Berry phase signatures 
are of the same magnitude as the nonadiabatic compo- 
nents. Berry phase effects in the AB-oscillations of the 
total conductance and the UCF's are hence masked by 
the regular, nonadiabatic contribution [23]. However, in 
this broad plateau region one still finds effects of the 
inhomogeneous magnetic field that can be ascribed to 
the nonadiabatic generalization of the Berry phase, the 
Aharonov-Anandan phase [24]. Our numerical results 
imply that to observe clear Berry phase effects such as 
the 'magic angles' found by Engel et al. [11] in the mag- 
neto conductance and UCF's of diffusive rings one has to 
go to the truly adiabatic regime given by condition (8). 

For a typical experimental AB setup based on cop- 
per rings with radius ro = 500 nm and £ = 15 nm [12] 
strict application of the criterion (8) corresponds to B- 
field strengths larger than 10 3 T. On the other hand, 
according to the condition Q <C l/y/(N), the opposite 
non-adiabatic regime ((Tn) ~ 1) is restricted to fields 
smaller than 0.1 T. In the broad intermediate B-field 
range, covering four orders of magnitude, the magneto 
conductance is expected to show at most signatures of 



the Aharonov-Anandan phase. 

Recently, imprints of Berry's phase in AB oscillations 
have been reported for holes in quasi-ballistic 2d GaAs 
rings with strong spin-orbit interaction [14] . A rough es- 
timate of the system parameters suggests that the exper- 
imental conditions may fulfill the adiabaticity criterion 
(8) with B replaced by an effective Rashba field strength. 

To summarize we studied spin-dependent quantum 
transport through 2d disordered geometries. We showed 
that the relevant parameter defining the adiabatic limit 
both for the total conductance and the quantum cor- 
rections scales with the square root of the number of 
scattering events. This can be cast into a generalized cri- 
terion for adiabaticity for both, ballistic and disordered 
systems. It appears as a severe obstacle for direct exper- 
imental observation of Berry phases in the conductance 
through diffusive metal rings. Our numerical findings 
indicate that elastic scattering due to non-magnetic im- 
purities in the presence of a spatially varying magnetic 
field exhibits features similar to those in systems with 
spin-flips associated with the scattering process as for 
magnetic impurities or Rashba spin orbit coupling. 
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